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We consider turbulent 4-wave interaction of two types of waves: acous- 

Q\ • tic waves (dispersion uj = k) and electromagnetic-type waves (dispersion 

q ■ S7 2 = m 2 + p 2 ). For large wave vectors (k 3> m), when the dispersion of 

d : 

CD ' EM-type waves becomes quasiacoustic, the stationary spectra are obtained 

CO 

following a standard Zakharov approach. In the small wave number region 
(k,p <C m) we derive nonlinear differential Kompaneets-type kinetic equation 
for arbitrary distributions of the interacting fields and find when this equation 

o 

CN ■ has Kolmogorov-type power law solutions. 

■ PACS 05.20.Dd, 05.30.Jp, 05.90. +m 

I. INTRODUCTION 

o 

In this work we consider 4-wave interaction of two types of waves: one with acoustic 



dispersion u = k, called massless field, and another with electromagnetic-type dispersion 
Q 2 = m 2 +p 2 , called massive field (u,Q are frequencies, k,p - wave vectors of the waves, 
the phase speed of the acoustic waves was set to unity). This is a fairly common case 
of wave interaction. It includes, for example, interaction of electrons with photons and 
electromagnetic fields with Langmuir waves. We assume that the interacting fields are 
scalar fields obeying quantum Bose-Einstein statistics. Then the collisional integral becomes 
a nonlinear functional in occupation numbers of both types of waves. We are interested in 
finding stationary solutions of such collision integral. In particular we will be looking for 
power law type solutions (n oc k a , N oc p 13 ) which can exist in the inertial ranges k, p ^> m 
or k, p <S m. 
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A standard procedure to find power law solutions of the kinetic equation is through 
Zakharov conformal transformations [[[]]. This method provides a simple rule for finding 
stationary power law solutions in the inertial region for differential scattering (when the 
change in energy or momentum in each scattering is small). For coupled kinetic equations 
the Zakharov transformation in the general cannot be done. Tractable interactions include 
(i) 3-wave interaction of waves with similar dispersion relations (using conformal transfor- 
mations), (ii) 3-wave interaction of high frequency and low frequencies waves [p]] (using a 
diffusion approximation) . 

In this work we consider coupled kinetic equations for 4-wave interaction of waves. In 
the "ultrarelativistic" regime, when the wave vectors of both particles is much larger than 
the mass k,p ^> m, both types of waves have similar acoustic-type dispersion and thus 
can be viewed as selfinteraction between quasiacoustic waves. The stationary spectra then 
can be found using conformal transformations. In the oposite limit of small wave vectors, 
k, p <C m, different dispersion laws of interacting waves break the homogeneity of the kernel 
of the kinetic equation and do not allow exact Kolmogorov-type solutions. In this case 
we expand the collision integral in small changes of energy in each collision. Then the 
collision integral may be simply written as a nonlinear diffusion equation. This method was 
first applied by Kompaneets || who derived a nonlinear equation for the interaction of the 
low frequency boson particles (photons) with nonrelativistic classical electrons which were 
assumed to be in thermodynamic equilibrium due to quick relaxation by the long range 
Coulomb forces. 

The "ultrarelativistic" regime, when both interacting waves have quasiacoustic dispersion 
requires a careful consideration - for exactly acoustic dispersion the weak turbulence theory 
may not be applicable (e.g., [0- |3J). Balk [JJ] and Newell and Aucoin || argued that in 
the multidimensional case the angular dispersion can justify the applicability of the weak 
turbulence theory to acoustic waves, while Zakharov and Sagdeev resorted to small 
corrections to the acoustic dispersion. We follow the latter approach: in the limit p ^> m 
the correction to the dispersion of EM-type waves will provide the dispersion required for the 
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final time of interaction of acoustic waves and will ensure that kinetic approach is applicable. 

The primary purpose of this work is to derive a low frequency diffusion approximation 
for the 4-wave interaction for arbitrary distribution functions of both interacting waves and 
to find Kolmogorov-type solutions in the inertial ranges. In case of quantum fields, the 
Kolmogorov-type power law solution may be possible only in the inertial ranges, where 
energies of the interacting particles are not close to the only quantity with dimension of 
energy - m and the occupation numbers satisfy n,iV > 1 or n,^ < 1. 

We start with a Hamiltonian of the system of two interacting waves in terms of wave 
amplitudes a k and A p 



H = J co k a k a* k d 3 k + J Vt p A p A* p d 3 p + 

7 / Tkpk>p>atA*a k >Ap>5(k + p - k' - p')d 3 kd 3 pd 3 k'd 3 p' 



(1) 



We assume the simplest form of the relativistically invariant matrix element (e.g., 0) 
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T k pk'p' 
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to ' k VtpU) ' k iVtpi 

Then the kinetic equations for occupation numbers n k , N p are 

\ J T'U>p>F[n, N]8{u k + tl p - u k , - Q^)S{k + p - k' - p') 

J 

where 
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d 3 k'd A p' (3) 



F[n, N] = N p ,n k , (n k + N p + 1) - n k N p {N p , + n k >) 



(4) 



II. APPLICABILITY OF THE KINETIC EQUATION 

For the applicability of the weak turbulence theory (based on random phase approxi- 
mation) we need that the typical interaction time ti nt ~ l/(TN to t), (where N to t = / d 3 pN p 
and T is a matrix element for amplitudes) be much larger than chaotization time for the 
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phases of waves t co m S i on = 1/Akv g (v g is the group velocity and the wave packet is centered 
on k ± Ak) and diffusion time t di ff = 1/ Ak 2 v' g . 

For a power law to ~ k a we have v' g ~ a(a — l)k a ~ 2 , so t co iu S i on ~ which is can be 

very short - of the order 1/uo. On the other hand the diffusion time t di ff = (Ak/k)~ 2 / (a(a — 
Which obviously unapplicable to acoustic turbulence. Still, if a ^ 1, then the diffusion 
time tdiff also can be very short. This case (nonacoustic turbulence) is applicable for small 
wave vectors: k <^m. 

^ = TN/kc ~ (5) 
t^nt m 

where we estimated T(k <C m) — g 2 /(km) , N = nk 3 . Thus the applicability of the kinetic 
equation in this case requires 

which for k < m reducesto n< 4. 
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On the other hand for semiacoustic waves with k ^> m we find uj = kc{ \ + we have 
= t d iff = (^) 2 ^ ~ ^2, then the applicability criterion is 

!ffl = ^ = M«, (7) 

t in i m 2 
(the matrix element in this limit is T = g 2 /k 2 ). This will eventually break at 

m 

Thus for applicability of a kinetic equation in a range m <C k <C /c maa; the same condition 
should be satisfied, ti€^. 

We can also estimate when the higher terms will be unimportant in the kinetic equation. 
The ratio of the second order terms of expansion in g 2 to the first order is typically 

[2]/[f]=sWAcW2^ (9) 
For acoustic turbulence k~^> m the frequencies are Ui ~ k, this gives 



[2]/[l] a = ^ 2 «lifn«^ (10) 

while for small wave vectors ui\ = m uj 2 = k 

[2]/[l] = g 4 n 2 k/m < 1 if n < ^ (11) 

Thus we conclude that for existence of inertial range where weak turbulence theory is 
applicable the occupation numbers should satisfy 

N p ,n k <&\ (12) 

III. INTERACTION OF PARTICLES AT LARGE WAVE VECTORS 

In the limit k,p ^> m ("relativistic particles") the dispersion law for the EM-type waves 
become almost acoustic 

The applicability of the weak turbulence theory to the interaction of acoustic wave in more 
than one dimension is a long standing question. In 1-D weak turbulence theory is not 
applicable since the interaction time between two waves is infinite, so that mutual influence of 
waves becomes large. In larger dimensions there were claims ( @, @) that angular dispersion 
may limit the interaction time sufficiently to allow weak turbulence approximation. Another 
approach is to allow for small corrections to the dispersion which will limit the interaction 
time of wave packets and allows us to implement methods of the weak turbulence theory 
regardless of angular dispersion [Q. 

Assuming the weak turbulence theory is applicable (Eq. (|12~1) ) we can find stationary 
solutions to the kinetic equation using Zakharov conformal transformations. In the limit 
k,p 3> m the dispersion laws uj oc k a and Q oc p a with a — 1, and the interaction coefficient 
in this limit is a homogeneous function of the order —2: 

T kpkV (\k, Ap, Ak', Ap') = A p T fc i2 3 (k, p, k', p'), p = -2 (14) 
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Zakharov transformations then gives stationary solutions 



Hp, Np oc LU 



5/3 



n Q , N Q (xu 



4/3 



(15) 



Here P and Q symbolize Kolmogorov spectra with constant flux of energy and constant flux 
of action. It is possible to check that collision integral converges for both spectra, so that 
the interaction is local. 

The sign of fluxes sign[P] > and sign[Q] < imply that two possible stationary spectra 
include the constant energy flux flowing to large k or constant flux of particles to small k. 
Realization of a particular spectra (Eqns flT5D depends on the boundary condition at the 
limits of the inertial range, e.i., by the location of the sinks and sources of energy and action. 
If the source of energy and of action are not at k = and k = oo correspondingly, then, 
depending on values of p and a, the integrals of energy or action converge or diverge. Using 
criteria for the stability of Kolmogorov spectra ||, we conclude that for sources located 
at k ,p 3> m it is Kolmogorov-in- action case that is realized (the action integral diverges 
at large k). In this case the turbulence is Kolmogorov in action - particles are flowing to 
small k (and may eventually be stored in Bose condensate). For k,p < ko,po a stationary 
Kolmogorov-in-action spectrum will form, while for for k,p > /co,Po there will be a selfsimilar 
"thermal wave" type solution carrying energy to large k and p: 



where it is eventually dissipated ||. Thus, in the "thermal wave" the mean frequency 
increases with time oc t^. 

This situation (flow of particles to small wave numbers, where stationary spectrum is 
formed, and flow of energy to large wave numbers) is somewhat similar to interaction of Lang- 
muir waves where in the low frequency region the dominant nonlinear process for Langmiur 
waves is the exchange of virtual ion sound waves which leads accumulation of particles with 
small k - Langmiur collapse. 




for k > h 



o 



(16) 



6 



IV. INTERACTION OF PARTICLES AT SMALL WAVE VECTORS 



A. Kompaneets approximation 

In the limit k,p <C m the change in energy in each collision is small, so we can expand 
the collision integral in small changes of the energies of particles: 



{tip' — {tip — A 



(17) 



where A to the first order in u/m is given by 

uj (uj (1 — n ■ n') + p • (n — n')) 
m + uj (1 — n • n') — p • n' 

Expansion of F[n, N] to the second order in A gives 

/ dn dN\ 

F[n, N] = A [N(N + 1)^ - n(n + 1)^- + 

\ OUJ OLO J 

ir +l fe +n(n+l) ^- 2 ^ (n ^ +i) 



(18) 



(19) 



In the following we concentrate on the kinetic equation for particles n. In the kinetic 
equations we next separate angular and momentum integrations 



I coU (k) = Ana / p 2 dp 
Jo 



^^A((*,p,„<)> 



where 



^ d2N » r /, r ^d 2 n dNdn, AT 

n{n + 1) 7w +N{N + 1) ¥^- 2 dn^ {N + n + 1) 



(20) 



Denoting the angle averages 
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fi\ 



IQir m 2 



r dtl p dfL n , 
h(p,k) = a J— — 



.... Is k (4k 2 -3km + p 2 ) 
A((fc,p,n') = — ^ — ^a 



h(p,k) = (7 J 



__A((*,p,n) = 



3 m 2 

,. 2 _ 2 A; 2 (2P + p 2 ) 



3 m 2 



(21) 



(22) 



We can rewrite collision integral in the form 



poo 

Icoii(k) = 4tcct / p 2 dp 
Jo 

1 / d 2 N 



(dN 

hM [m n{n + 1) ~d^ N{N + Vj 



+ N{N+1) 



d 2 n 



dn 
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d 2 ui dil duj 



(N + n + 1) 



(23) 



There is a also a constraint condition on the equation ( |2"3"D - it sould conserve total 
number of particles, which is an integral of motion of the Hamiltonian (Eq. [I]): 



J d 3 kI coll (k) = 



(24) 



Partially integrating and collecting terms with different powers of distribution function we 
find 



dudQNn 
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(25) 



where 



^l 



2 dk 2 dp . 2 dk 2 dp 

k ~X P I\i l2 = k Z —p Z —I 2 



duj dVt 



du dVt 



(26) 



Each term in Eq. 
%\ and i 2 : 



5|) should be equal to zero. This gives a relation between the coefficients 



1 / di 2 di 2 

2 vM ~ 



(27) 



This is an important relation that is one of the main results of the work. It relates the coef- 
ficients in the diffusion expansion of the collisional integral independently of the particular 
form of the interaction. 

Using relation fl2"7j) we can simplify considerably the collision integral (Eq. |3[). Partially 
integrating we find that collision integrals take a particularly simple diffusion type form: 



S 



I coll (k) = 2na J P 2 dp-^ {^h [n'N(N + 1) - N'n(n + 1)] 

/O 



QI 2 [n'N(N + 1) - N'n(n + 1)] 



(28) 



For massless particles n and massive particles N correspondingly. 

Equations (|28|) are the diffusion-type equations for the low frequency interaction of two 
types of quantum waves - "massless" acoustic waves "n" and "massive" EM waves N. To 
reduce then to the more familiar Focker-Plank form we introduce coefficients 



A x = Ait jp 2 dpN(N + 1), A 2 = 4ir J k 2 dkk 2 n(n + 1) 



J p 2 d P N(N+l)p 2 
j p 2 dpN(N+l) 

J p 2 dpN'p 2 
' j p 2 dpN(N+l)p 2 



Q 
B 2 



Jk 2 dk^n(n + 1) 
J k 2 dkk 2 n(n + 1) 
/ k 2 dkk 2 n' 



J k 2 dkk 2 n(n + 1) 

Jp 2 dpN' P 2 Jp 2 d P N(N+i) _ r _ I k 2 dk^n! J k 2 dkk 2 n(n + 1) 
jp 2 d P N' jp 2 dpN(N+i)p 2 ~ ' 2 ~ JkMkk 2 n' J k 2 dk^n(n + 1) 



(29) 



Note, that in the thermodynamic equilibrium N' = —N(N + 1) /T, so that Bi, B 2 = 1/T 
(T is equilibrium temperature) and Ci,C 2 = 0. Coefficients v4 lj2 and u, Q are complicated 
even for thermodynamic equilibrium. In the classical limit they give A\$ = N to t, n to t - total 
number of particles, uo 2 = 3mT -average < k 2 > and flo — 5T/4. 

Using notations (|29|) the collision integrals may be written as 



-fcollM 



O"0 



Icoll(fi) = y^2 



VIqC 2 



+ fi ) ( ^ + B 2 N{N + 1) (l + 



(30) 



B. Stationary solutions 

To find stationary power law solutions to the nonlinear equations (^) we assume that 



Ci,2 ~ (in this case the resulting solution are sums of power laws (see Eqns. (|35|-|3"6|)). 
In case of C\ the dominant contribution to the number density comes from the short wave 
vector limit of the thermodynamic equilibrium (see Eq. (|36[) ), for which we can assume 
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G\ ~ 0. The coefficient C 2 on the other hand is generally nonzero, but if the maximal 
energy of particle n (cutoff of the power law) is much less than the mass m then the term 
with C 2 is approximately p m ax/m - much smaller than unity. 

Under assumption C\^ 2 ~ the Kompaneets-type equations for the occupation numbers 
n and N become 



dn cr ^ 1 d 



dt 2 2 Jttdn 



Q(Q + Q ) (^L + B 2 N(N + 1) 



(31) 



and the the stationary equations || — — 



dt 



should satisfy 



u 4 {u 2 + ul) (ri + B in (n + 1)) = P 1 
Q(Q + ft ) (N' + B 2 N(N + 1)) = P 2 



(32) 

The physical meaning of the constants Pi and P 2 is the total particle flux in the wave vector 
space. 

Eqns (B^) are nonlinear Riccatti-type differential equations. Using a conventional sub- 
stitution N'(x) — > U'/U they can be reduce to linear equations of the form 



U" + U'- G{x)U = 



(33) 



It is not generally integrable in quadratures, but in some limits it may have a solution which 

is a sum of power laws ||. [] The best way to search for possible power law solutions is to 

-1 



use the substitution U(x) — > W exp{yj yG(x)J dx} which takes a correct account of the 
leading power law term. 

Power law solutions of the equation ([32] ) are possible in the inertial regions - where there 
is no quantity with a dimension of length: these are regions u much less or much larger than 
ujq, Q much less or much larger than Qq and N, n much less or much larger than 1. 



For the equilibrium distribution of heavy particle N various approximations for n in the limit 
u> S> u>o an d N <C 1 are summarized by |f[0(| . 
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In the linear regime Af,n < 1 the stationary solutions are 

n(X i 3? f 1 , — 2^ duJ 

J LU d (iU 2 + UJq) 

N oc / --^ rdQ 34 

In the strongly nonlinear (N, n ^> 1) limit the leading orders in power-laws are 



n,„ oc < 



-^T + -W- for Co> <C UJn 

ujqw^ Bu ^ u (35) 

Pi , 3 



+ 2^ for ^ » ^0 



)7?o T7^T7T fo r ^ < ^0 
(36) 
T^T - f#4 for O » 

Given the solutions (B6|), we can also calculate the coefficients Qq and ujq\ 







8m 



for n oc w 2 



— for n oc uj 



12m 

U0=Pmax/V3 (37) 

where k max and p maa ; are the upper cut-offs limits for the corresponding power laws. 

We can verify that the integral in the expressions for the coefficients (Eq. (|29|) converge 
at k,p — 0. Note also, that for massive particles N, the leading term in the small wave 
vector limit corresponds to the low frequency asymptotics of the thermodynamic equilibrium 
spectrum. 



Equations ( j35| - |3~7l) give the selfconsistent Kolmogorov-type solutions to turbulent inter- 
action of two types of waves. They are qualitatively shown in Fig. (0). In all cases the 
particles are flowing to small k,p with a constant flux in phase space. If the sources of 
particles are located at k,p ^> m (Fig. |2] (a)), so that ujq « m, then the low energy tail of 
the distribution of massless particles (u <C u ) is n oc uo~ 2 . Alternatively, if the sources of 
particles are located at k,p <C m (Fig. || (b)), then the low energy tail of the distribution 
of massless particles is n oc 
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V. CONCLUSION 



We have considered turbulent 4-wave interaction of two types of waves with different 
dispersion laws: acoustic and electromagnetic- type. Both types of waves were assumed to 
obey a quantum boson statistics. We found possible stationary Kolmogorov-type spectra in 
the inertial ranges k, p 3> m and k, p <C m. To reach the stationary solutions the particles 
should be generated at finite momenta k and p and flow to small wave numbers. In the 
"ultrarelativistic limit" k, p 3> m, the turbulent interaction of two types of waves resembles 
self interaction of semiacoustic waves with the dispersive correction ensuring the applicability 
of the weak turbulence limit. In the small wave vector limit k, p Cm the collisional integral 
may be reduced to diffusion- type differential equation. In this limit, the small wave number 
asymptotic of massive particles is dominated by the thermal equilibrium tail N oc 
while the small wave number asymptotic of massless particles is steeper: n oc oj~ 2 -t- u; -3 . 
Arguably, the massless particles n will eventually reach k = to be absorbed in the Bose 
condensate.0 

I would like to thank Alexander Balk for useful comments. 
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FIGURE CAPTIONS 

FIG. I. Stationary power law solutions of the system (0). (a) source of particles is 
located at k,p ^> m, (b) source of particles is located at k,p <^ m 
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